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A theory of low temperature (T) electron spin resonance 
(ESR) in half-integer spin antiferromagnetic chains is devel- 
oped using field theory methods and avoiding previous ap- 
proximations. It is compared to experiments on Cu benzoate. 
Power laws are predicted for the line-width broadening due 
to various types of anisotropy. At T — > 0, zero width ab- 
sorption peaks occur in some cases. The second ESR peak in 
Cu benzoate, observed at T < .l&K, is argued not to indi- 
cate Neel order as previously claimed, but to correspond to a 
sine-Gordon "breather" excitation. 

PACS number 75.10.Jm 



In ESR a static magnetic field is applied and the ab- 
sorption of radiation polarized along the field direction 
is measured, as a function of frequency. The absorption 
intensity is proportional to the Fourier transform of the 
correlation function of the (zero wave-vector) total spin 
operator, G[t) =< S^{t)S~{{)) >, when the field is along 
the z-axis. (Here S = ^ Sj.) We write the Hamiltonian: 



Ti. ~ Ti.Q + Ti' + Tiz 



(1) 



where Hq is the isotropic Heisenberg Hamiltonian, 
J^Si • Si+i, Hz is the Zeeman term, ~H^S^, and 
Ti' represents various possible small anisotropic terms. 
(We set g^lB = !•) When W = 0, [n,S-] = HS-, 
implying that the absorption spectrum consists only of 
a ^-function peak at the Zeeman energy, H. The shift 
and non-zero width of this peak are caused by the small 
anisotropic terms in the Hamiltonian, Ti.' . 

In this paper we develop a new approach to calculating 
the ESR intensity for one dimensional antiferromagnets 
(ID AF's) of 1/2-integer spin. Using bosonization and 
the standard Feynman-Dyson self-energy formalism we 
are able to avoid making assumptions about the relax- 
ation function as in the previous treatment 

We then calculate the width and shift perturbatively, 
avoiding previous Hartree-Fock approximations which 
are generally invalid in ID. This perturbation theory gen- 
erally breaks down due to infrared divergences at low T 
where universal scaling functions give the shift and width, 
whose calculation requires more powerful methods, based 
for example on exact integrability. We study the case 
where H' corresponds to a staggered field We show 
that the width is proportional to {H/TY and the shift to 



[H /TY at intermediate temperatures. The predicted de- 
pendence on T, H and field direction is shown to agree 
with experiments on Cu benzoate. We argue that the 
low temperature ESR experiments are observing excita- 
tion of a sine-Gordon breather above the groundstate, 
rather than Neel order as previously claimed. 

We bosonize the spin operators, with Tio corresponding 
to a free boson Lagrangian, (1/2) J dxdfj^cpd'^cj) and: 



S- 



/2Trdx 



+ constant (—I)-' cos(v 27r(/)) 



iV2TT<1 



[constant cos V27r(/) + C(— l)-*]. 



(2) 



Here, the fields may be decomposed into left and right 
movers a.s (j) — (pL + (jfR and 4> = (pL — 4>r and we 
set the spin- wave velocity to 1. Noting that Tiz = 
— {H/\/2T:)d4>/dx, we see that Hz may be eliminated 
by the field redefinition. 



(3) 



{H/^/2-k)x 



This shift must be applied to the bosonization formulas 
of Eq. (|^). Upon Fourier transforming the spin opera- 
tors, this means that some momenta get shifted by ±iJ. 
Writing the low-momentum parts of the spin operators 
in terms of left and right spin currents, (Jl -I- Ji?), we 
see that the ESR absorption intensity can be written in 
terms of the Green's functions of these operators. The 
effect of shifting (j) is to shift the Fourier modes of the 
currents as j|(fc) ^ J^{k + H), J^{k) ^ J^{k- H). 
Note that the bosonized version of Ho+Hz is apparently 
independent of H and hence SU(2) invariant; it is only 
the mapping from the lattice spin operators to the field 
theory which depends on H . (Actually, it is known that 
the correlation exponents vary with H, corresponding to 
a change in the compactification radius of (f>. However, 
this can be seen to be irrelevant to ESR, which probes the 
finite energy H . The effective SU(2) symmetry is man- 
ifested by the zero line- width of the ESR peak.) The ESR 
intensity is determined by the Green's functions oi j^, 
which are given by exponentials of 4'l,r- On the other 
hand, the Green's functions of J£ ^ are easier to deal 
with since these fields are linear in 4'l,r- In particular, 
we may use the standard Dyson result to write the re- 
tarded Green's function for (p in terms of the self-energy, 



n(g,^,T): G 



/c~ — n] . This is a very useful 
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formula for ESR because, in the limit where we may treat 
TC perturbatively, the shift and width are simply given 
by the real and imaginary parts of 11 (iJ, H). Using the 
effective SU(2) symmetry of the bosonized theory we may 
express the ESR intensity in terms of the Green's func- 
tions of J£ ^. Note that the validity of Dyson's formula 
is a highly non-trivial result of the structure of pertur- 
bation theory to all orders, resulting from the multiple 
insertions of one-particle irreducible diagrams and free 
propagators. In our new approach, this essentially re- 
places the previous assumptions ^ about the relaxation 
functions. 

We consider the particular example of a transverse 
staggered field: 



(4) 



This term arises, with h H, from either a stag- 
gered off-diagonal component of the g-tensor or from 
a Dzyaloshinskii-Moriya (DM) interaction |2|; both of 
these occur in Cu benzoate. The bosonized interaction 
is Ti.' = hC / da; cos (\/27r0), which is unaffected by the 
field redefinition. In the case H = 0, TCq + Ti.' is invari- 
ant under rotation about the z-axis by tt. It then follows 
that the Green's function giving the ESR intensity can 
be expressed as G+- = G=""' + G^y, since G^'' = 0. (Here 
G"'' = -i < [S^'^S^] >.) We now use the SU(2) sym- 
metry of Ho to prove that G^^ for a staggered field in 
the x-direction is the same as G^^ for a staggered field 
in the z-direction and G^^ for a staggered field in the x- 
direction is the same as G^^ for a staggered field in the 
y-direction. As argued above, this SU(2) symmetry re- 
mains present in the bosonized version of Tio + Ti-z for 
non-zero H. Thus we have succeeded in expressing the 
ESR intensity in terms of < (j)(f) > although we appar- 
ently must consider 2 different forms of H' . In fact these 
2 different H's only differ by the interchange of (p ^ <j>, 
corresponding to ^Jr- Thus, the contribution to 

the ESR intensity, cx< J£J£ > -I- < J^J^ >, is identical 
in both cases. Thus our formula for the ESR intensity, 
I[ijj,H,T) becomes: 



/ cx —uj\m 



i?2 -n(w,iJ,T) 



(5) 



This fundamental formula replaces the more ad hoc 
ones generally used in ESR theory. It will be partic- 
ularly useful when H.' can be treated as a small per- 
turbation, resulting in a small value of Il{H, H) / H <C 

T . In this case Eq. (|^) predicts an approximately 
Lorentzian line-shape with shift JieI{{H,H)/2H and 
width -Imn(7J, H)/2H. 

Note that /i is a relevant coupling constant of scal- 
ing dimension 3/2. It produces an excitation gap, A cx 
ji/3/j2/3^ Thus perturbation theory in h will only be 
valid at high T; infrared divergences occur for T < A. It 



follows from general scaling arguments that the shift and 
width both have the form Tf,{A/T, H/T), for 2 different 
scaling functions fi. We first consider the case T 3> A 
where we may use lowest order perturbation theory in h; 
we consider the low T case later. The first non- vanishing 
term is 0{h?). We must amputate the external lines from 
the Green's function < (/)exp[i'y27r0] exp[— z\/27r(/)]0 >. 
This can easily be done by Taylor expanding the expo- 
nentials and gives two terms corresponding to Feynman 
diagrams in which the external lines are attached to the 
same or different vertices. This gives: 



Il{H,H) « 2Tr{Chf[G\H, iJ) - G'(0, 0)], 



(6) 



where G' represents the retarded Green's function of 
sin(-\/27r0) in the free theory. We also performed a field- 
theory calculation using the traditional approach j^] to 
find the same result as above, in the lowest order of per- 
turbation. Of course the higher order perturbative terms 
from our new formulation, and in particular the non- 
perturbative low T behavior, go beyond the traditional 
approach. 

G' has been evaluated by Schulz In general for 
a normalized primary field with left and right scaling 
dimensions xl — xr = x/2 the retarded Green's function 
is given by, using u!± — lu ±q, 



G(uJ,qU 



(7) 



F^{e) EE r(a;/2 - ie/47r)r(l - 2:)/r(l - x/2 - ie/An) (8) 

where F is Euler's Gamma function. For the staggered 
field 1/2. Also considering the limit i? -C T", we 

may Taylor expand -^^^(e). Although the proportionality 
factor, G, in the bosonization formula, Eq. (|^) is not 
universal it was recently evaluated exactly using Bethe 
ansatz results [Q. Due to a marginally irrelevant opera- 
tor in the Hamiltonian, so far ignored, G^ is effectively 
proportional to InT (for T ^ H). This finally gives the 
results for the shift, Suj and width rj: 



Slu : 

V ■■ 



A2596{\n J/T)Jh^H/T^ 
.68570l(\n J /T)Jh^/T^. 



(9) 



The effective staggered field, h, is determined by a linear 
combination of the staggered component of the g-tensor 
and the DM interaction. It takes the form h = cH where 
the proportionality constant, c, is a strong function of 
the direction of the applied field. Thus the shift and 
width should scale as essentially (H/T)^ and (H/T)'^ re- 
spectively, with strong field- direct ion dependence, in the 
temperature range A, <§C T <§C J. (In experiments, 
usually the frequency is kept fixed; the frequency de- 
pendence corresponds to the i7-dependence.) Note that 
the width (and shift) increase with decreasing T. The 
increase of the shift due to the staggered g-tensor was 
already discussed by Nagata jsj, although our result is 
quite different from his. 
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FIG. 1. The field direction dependence of the ESR 
hnewidth in the ac-plane at frequency 35GHz ^ compared 
to Eq. (|). 
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FIG. 2. The temperature and frequency dependence of the 
ESR linewidth for H|| c" i], compared to Eq. (jgl). 



Such a peculiar behavior of the width was actually ob- 
served in Cu Benzoate We choose the DM vector, 
{Da" , Dc") = (.13, .02) J, which seems consistent with 
most of the experimental data on Cu benzoate Q . Fig- 
ure 1 and 2 shows the dependence of the width on field 
direction, T and H. (There we replaced the logarith- 
mic factor in eq. (^) by unity.) We see from Fig. 1 
that while the direction dependence and magnitude of 
the width are roughly reproduced there appears to be an 
additional smaller contribution to the width with weaker 
dependence on field direction. We expect such a contri- 
bution from exchange anisotropy (discussed below); there 
may be other contributions as well. In Fig. 2, the tem- 
perature dependence of the linewidth for H||c" is shown 
for several resonance frequencies (ie. field magnitude H) . 
We see that the dependence of the width on H and T are 



well reproduced by our theory on the staggered field ef- 
fect, although an additional contribution to the width, 
depending less strongly on H and T, is evident. 

Another possibility is that the increasing width with 
decreasing T is related to the onset of Neel order at low T 
due to interchain couplings ignored in our theory. How- 
ever, there is no evidence for Neel order from neutron 
scattering or susceptibility measurements ^ . (As we ar- 
gue below, the interpretation of ESR at still lower T as 
evidence for Neel order seems incorrect.) Furthermore, it 
seems unlikely that the peculiar dependence of the width 
on the three different variables, T, H and direction, cap- 
tured by our purely ID theory, could be explained in this 
way. 

For T < A, the above perturbative analysis breaks 
down. In the extreme low T limit, T ^ A wc can 
nonetheless make some statements about the ESR in- 
tensity based on general principles and on exact results 
on the sine-Gordon model The simple picture of a 
Lorentzian line shape characterized by just a shift and 
a width is no longer very useful. Instead at very low T, 
ESR experiments should, in principle, be able to resolve 
some of same the excitations above the groundstate as 
seen in neutron scattering. [Related observations were 
made independently by Boucher et al. recently [|To|.] 
Using the SU(2) rotation trick discussed above, we see 
that the radiation can create the excitations of the sine- 
Gordon model produced from the groundstate by the 
scalar field, (p. This includes the first breather as well 
as a multi-particle continuum. Thus the ESR intensity, 
/(w) at T = should contain a zero width peak at the 
energy of the breather with momentum H, 



v/i/2+A(/l)2, 



(10) 



in addition to a continuum at higher energies. Calcula- 
tions of the (j) form factor |^ based on the exact integra- 
bility of the sine-Gordon model indicate that the inten- 
sity of the multi-particle contribution is very small and 
that the spectrum is dominated by the first breather peak 
at zero temperature. This breather peak can only get 
broadened by collisions with thermally excited particles 
so we expect that its width should vanish as exp[— A/T]. 
As the temperature is raised other contributions to the 
ESR intensity should appear corresponding to additional 
processes involving thermally excited particles. 

Experimentally a noticeable change in ESR spec- 
trum occurs at T of 0{A). In particular a two-peak 
structure is observed for an intermediate range of T. At 
lower T only the higher energy peak survives. A pos- 
sible interpretation of this behavior is that the higher 
energy peak represents excitations from the groundstate, 
perhaps primarily the first breather whereas the lower 
energy peak, which disappears with decreasing T, cor- 
responds to other processes involving thermally excited 
particles. This hypothesis can be tested by fitting the 
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shift to the breather energy of Eq. (|^. Note that this 
leads to a characteristic direction dependence of the shift 
due to the strong direction dependence of the effective 
staggered field, h and hence of the gap, A. Using the 
previously determined parameters, we obtain the depen- 
dence of resonance field (i.e. shift) on direction which 
is compared to experimental data at low T in Fig. 3. 
The agreement looks quite good. On the other hand, 
the observed width of this low T peak appears to go to 
a non-zero value at low T, which appears inconsistent 
with the breather interpretation. This may be due to an 
additional broadening mechanism which is effective at 
low T such as quenched disorder. Further experiments 
at higher field would help to clarify the situation. In 
any event, the previous interpretation of the low T peak 
as a signal of Neel order seems unjustified. Instead it 
can perhaps be explained entirely within the ID theory 
taking into account the non-trivial evolution of the ESR 
intensity with T and H . This would then resolve the con- 
tradiction between the claimed observation of Neel order 
from ESR and its non-observation in neutron scattering. 
It is worth emphasizing that the staggered field, by pro- 
ducing a gap and a staggered moment, tends to suppress 
the interchain coupling effects. 
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FIG. 3. The resonance field at very low temperature 
for various field directions in ac-plane. The experimental 
data were taken from the "Antiferromagnetic Resonance" in 
Ref. at 0.41K, and the theory refers to the lowest breather 
excitation at zero temperature. 



from scaling that rj cx Tf{5,H/T), where |M| cx 5 and 
/ is a scaling function. Using Eqs. (0), (||) in the hmit 
of a; ^ 2, and including the effect of the large isotropic 
marginal operator gives the width: 



7?(x(<5/Jf(in(j/r))2r. 



(11) 



Such a linear T-dependence of the width is observed ap- 
proximately, over an intermediate range of T, in a variety 
of quasi-lD antiferromagnets As far as we know, 

ours is the first derivation of this behavior from first prin- 
ciples. 
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We have also analyzed the contribution to the width 
from exchange anisotropy (and/or dipolar interaction) 

'12 j SfMaf:iSj_^i using field theory, based on our new 
formulation as well as the traditional one ||l|. We em- 
phasize that it is not necessary to use a Hartree-Fock 
approximation to calculate the Green's function, as has 
been commonly done in the past, even when the opera- 
tor is quadratic in the lattice spin operators. Indeed such 
Hartree-Fock approximations are generally qualitatively 
wrong in the ID case. Since H' is marginal it follows 
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